The dipole picture in DIS: saturation and heavy quarks 
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^«— ^ ' We discuss the description of the proton structure function within the dipole factorisa- 

tion framework. We parametrise the forward dipole amplitude to account for saturation 
as predicted by the small- 2: QCD evolution equations. Contrarily to previous models, 
J5 ' the saturation scale does not decrease when taking heavy quarks into account. We show 

^ , that the same dipole amplitude also allows to reproduce diffractive data and exclusive 

vector meson production. 
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In these proceedings [T] we shall concentrate on Deep Inelastic Scattering (DIS) at small 
X. In this regime, the photon-proton cross-section can be factorised as a convolution between 
the wavefunction for a virtual photon to fluctuate into a quark-antiquark pair and the 
interaction T between this colourless dipole and the proton: 

aL.T{x,Q^)=27rRlJ2 I d?rdz |^L,T(r, z; g')|' T(r,x), (1) 



where the factor iirBl arises from integration over the impact parameter. The photon wave- 
function can be computed from perturbative QED and we are left with the parametrisation 
f~^ ' of the hadronic dipole amplitude. To that aim, we usually rely on the observation that the 

^SJ , small-x DIS data satisfy geometric scaling [2], meaning that, instead of being a function of 

^D ' both Q^ andx, they appear to be a function of r = \og{Q'^ / Ql{x)) = logi^Q"^ /Ql) — \\og{l/x) 

^— ; , only, where Qs{x) is known as the saturation scale. Since r ^ 1/Q in ([T]), this property sug- 

C^^ ' gests that the dipole amplitude is a function oirQs{x) only. 

^^ , Since the small- a; domain extends down to small Q^, the dipole amplitude is sensitive 

to the unitarity bound T < 1. There are two broad classes of models which differ by their 
way to implement that boundary. The first approach is to use an eikonal form as initially 
proposed in [3|, followed by more precise analysis to incorporate DGLAP evolution and 
masses for the heavy quarks [4]. 
;J] ' The second approach, that we follow through these proceedings, is to use predictions 

5^ , directly from the Balitsky-Kovchegov equation describing the QCD evolution to small x. It 

resums the BFKL logarithms of \/x and satisfies unitarity by including saturation effects. 
In contrast with the eikonal models which include it by hand, it has been proven [5] that 
the solutions of the BK equation satisfy the property of geometric scaling. More precisely. 
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In this expression, 7c, A and k are obtained directly from the BFKL kernel. The first term 
in the exponential, surviving at asymptotically small x, satisfies geometric scaling, while the 
second term violates geometric scaling and describes how it is approached when x decreases. 
Geometric scaling is thus respected when the second term can be neglected i.e. in a window 
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model 


7c 


Vc 


xo 


Rp 


x'/n 


IIM 


0.6275 


0.253 


2.67 10-^ 


3.250 


wO.9 


IIM+c,b 


0.6275 


0.195 


6.42 10-' 


3.654 


1.109 


new fit [6] 


0.7065 


0.222 


1.19 10-^ 


3.299 


0.963 



Table 1: Values of the parameters and x^ per data point 
for (i) the original IIM model, the IIM model with heavy 
quarks and fixed 7c and (iii) the new, adapted, model. 



extending up to z = ^KAlog(l/x), which extends beyond the saturation momentum itself. 
This is an important message one learns from the BK equation: the saturation effects are 
not only important for rQs > 1 (or Q^ < Q^); they do affect the physical amplitude at 
larger values of Q^, in a window growing like •\/log(l/x) and where the dipole amplitude T 
is significantly smaller than 1. 

We thus can use ([2]) to 
parametrise T at small dipole 
sizes and match it continu- 
ously with an expression of 
the form 1 — exp(— (az + b)^), 
describing the solutions of the 
BK equation in the deep sat- 
uration domain. The param- 
eters A, Xo and Rp are then 
fittectl to reproduce the latest 
HERA measurements of the 

inclusive proton structure function for x < 0.01. This method has been successfully ap- 
plied by lancu, Itakura and Munier (IIM) [7], as shown in the first line of Table [TJ where 
the sum over quark flavours in H]) only account for three light quarks. 

One of the general issues of these mod- 
els is that, in both classes of models, once 
the mass of the heavy quarks is taken into 
account in ([T]), the saturation scale drops 
down by a factor « 2 (« 500 MeV instead 
of w 1 GeV). This is illustrated by the sec- 
ond line of Table [U where we see that, once 
the contribution from heavy quarks is in- 
cluded, the quality of the fit becomes poor 
and xq decreases severely. 

Recently [6], I have shown that it was 
possible to accommodate, for the first time, 
the IIM model to include heavy quark con- 
tributions without having the inconvenient 
that the saturation scale goes down. The 
underlying idea is to allow the slope jc to 
become a free parameter of the fit. As 
shown on the third line of Table [U this does 
not only brings the parameters closer to the 
original IIM model, especially xq as we will 
comment further later on, but it also results 
in a much better x^ ■ To obtain the parame- 
ters mentioned in TablefTJ one has restricted 
the Q2 range to Q^ < 150 GeV^, though the 
parameters are stable when we vary this limit. Note however that we do expect corrections 
from resummation of the DGLAP logarithms at larger Q^ . 




Figure 1: Description of the Fj' HERA data. 



^7c and k are fixed to the value predicted from the leading-order BFKL kernel. 
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limit of geometric 
soaiing window 



Figure [T] shows how well the F^ data 
[5] are reproduced. More interestingly, 
we have compared on Fig. [2] the satura- 
tion scale obtained in this new paranietri- 
sation with the one of the IIM model 
(bottom lines). They are clearly of the 
same order, showing that it is possible to 
include the heavy quarks in the dipole 
picture and at the same time keep a 
dipole amplitude with a saturation scale 
around 1 GeV at HERA. We also observe 
on Fig. [21 (upper lines) that the data 
up to Q^ = 5 — 7 GeV^ he inside the 
geometric-scaling window and arc thus 
sensitive to the physics of saturation. 

Now that heavy quarks are properly 
integrated into the picture, we can have 
predictions for the charm and bottom structure functions. One sees from Figs. [3] that we 
achieve a good description of those data. This figure also shows the predictions for the 
longitudinal structure function where the model once again agrees with the data. 



Figure 2: Saturation line and extension of the 
geometric scaling window with HERA data. 
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Figure 3: The charm structure function (left plot), bottom structure function (top right) 
and longitudinal structure function (bottom right). 

With this new model for the forward dipole amplitude that includes heavy quarks and 
remains compatible with small-x QCD evolution, we can also look at exclusive processes. As 
we show now, one of the major power of the dipole picture is that it allows, with the same 
parametrisation of the dipole amplitude, to describe both inclusive and diffractive processes. 

The first of those measurements we shall consider is the diffractive structure function. 
Generally speaking, this quantity probes correlations inside the proton and is sensitive to 
the square of the dipole amplitude. However, considering only fluctuations of the photon 
into a colourless qq state only allows to describe the limit of large /3 (or the limit small 
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diffractive mass Mx <C Q). To go to smaller values of /? (keeping Xpom = x/P <^ 1), we 
also have to consider the radiation of one additional ghion from the initial quark-antiquark 
pair. We can then show that the interaction between the resulting qqg pair and the proton 
can also be expressed in terms of the dipole amplitude T, leading again to a contribution 
proportional to T^. Based on the new parametrisation [^, it has been shown [9] that the 
HERA measurements of the diffractive structure function are well reproduced. 

The final observable we shall consider is the pro- 
duction of exclusive vector mesons or on-shcU photons 
(DVCS). The differential cross-section da'^'P^^P/dt can 
also be factorised, this time as a convolution between a 
vertex for the virtual photon to fluctuate into a qq pair, 
the interaction between that dipole and the proton, and 
the vector-meson wavefunction to account for the final 
state. The main difference with ([T]), beside the pres- 
ence of the vector-meson wavefunction, is that one has 
to account for the momentum transfer dependence of the 
dipole scattering amplitude. An intuitive method is to 
make a Fourier transform and go to impact-parameter 
space. We shall rather use the result of studies of the 
full BK equation including its momentum-transfer depen- 
dence. It is predicted [10] that the saturation scale is 
constant at small t and increases like |t| at large t. Intro- 
ducing one parameter to implement that dependence and 
a second one to describe the proton form factor (taken 
of the form exp(— 6|t|)), we have reached [TTj a successful 
description of the differential and total cross-sections for 
exclusive productions of p, cj) and J/^ mesons, as well as the for DVCS measurements. The 
description of the p-meson differential cross-section is given as an example on Fig. [51 
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Figure 4: Description of the p- 
meson differential cross-section. 
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